MEMORANDUM 

BM588S-PR 

FEBRUARY  1S69 


INVERSE  FUNCTIONS  OF  THE 
PRODUCTS  OF  TWO  BESSEL  FUNCTIONS 

William  Sollfrey 


PREPARED  FOR: 

UNITED  STATES  AIR  FORCE  PROJECT  RAND 


PE3T  AVAILABLE  COPY 


MEMORANDUM 

RM-5886-PR 

FEBRUARY  1969 


INVERSE  FUNCTIONS  OF  THE 
PRODUCTS  OF  TWO  BESSEL  FUNCTIONS 

William  Sollfrey 


This  research  is  supported  by  the  United  States  Air  Force  under  Project  RAND— Con¬ 
tract  No.  F41620-67-C-0045— monitored  by  the  Directorate  of  Operational  Requirements 
and  Development  Plans,  Deputy  Chief  of  Staff,  Research  and  Development,  Hq  USAF. 
Views  or  conclusions  contained  in  this  study  should  not  he  interpreted  as  representing 
the  official  opinion  or  policy  of  the  United  States  Air  Force. 

DISTRIBUTION  STATEMENT 

This  document  has  been  approved  for  public  release  and  sale;  its  distribution  is  unlimited. 


■7 it 


RAPID 


1700  MAIN  »T  •  I A  N  f  A  MONICA 


CAt  IIOINII 


This  study  is  presented  as  a  competent  treatment  of  the  subject,  worthy  of  pub¬ 
lication.  The  Rand  Corporation  vouches  for  the  quality  of  the  research,  without 
necessarily  endorsing  the  opinions  and  conclusions  of  the  authors. 


Published  by  The  RAND  Corporation 


iii 


PREFACE  AND  SUMMARY 

This  Memorandum  is  a  part  of  RAND's  continuing  study  on  electro¬ 
magnetic  propagation  as  affected  by  atmospheric  scattering.  Atmos¬ 
pheric  properties  have  been  inferred  from  scattering  coefficients, 
whose  calculations  involve  finding  certain  inverse  operators.  Special 
cases  of  the  inverse  function  of  the  product  of  two  spherical  Bessel 
functions  have  been  found  recently  by  other  writers  as  ^2  hyper¬ 
geometric  functions.  This  Memorandum  gives  the  general  expression  as 
the  derivative  of  a  product  of  spherical  Bessel  functions.  These 
results  may  have  application  in  various  scattering  configurations, 
and  also  to  the  classical  physics  problem  of  determining  the  potential 
from  the  phase  shift. 
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INVERSE  FUNCTIONS  OF  THE  PRODUCTS  OF  TWO  BESSEL  FUNCTIONS 

A  recent  article  considers  the  problem  of  finding  the  inverse 
function  of  the  product  of  two  spherical  Bessel  functions.  The  results 
are  given  in  terms  of  ^F2  hypergeometric  functions.  The  inverse  func¬ 
tion  may  be  defined  by  the  equation 

I  (1) 

0 


Then  Ref.  1  gives  formulas  in  terms  of  ^F,,  functions  for  m  =*  0,1,2, 

and  explicit  results  for  g^  Q,  g1  Q,  82jo*  80,1*  and  g0,2*  The  authors 

conjecture  that  the  higher  order  inverse  functions  can  be  found  by 

the  methods  presented  in  Ref.  1. 

This  problem  has  been  treated  in  the  classical  literature.  The 

inverse  function  g^  q  was  found,  in  a  form  differing  by  an  integration 

by  parts,  by  Bateman, ^  and  the  higher  order  forms  by  Fox. ^  All 

(4) 

these  results  are  available  in  a  standard  text.  The  general  and 
quite  simple  formula  is: 


g,  _00  = 
l  ,m 


8x' 

n 


dx 


niwWx) 


(2) 


where  n^(x)  denotes  the  spherical  Neumann  function.  Equation  (2)  re¬ 
duces  to  all  the  special  cases  treated  in  Ref.  1  except  for  g.  -, 
for  which  it  differs  by  a  constant.  Since  it  may  be  easily  verified 
that  and  j^2  are  themselves  orthogonal  over  the  range,  the  in¬ 
verse  function  is  not  unique  to  an  arbitrary  additive  constant. 
Furthermore,  if  the  spherical  Bessel  functions  are  replaced  by  the 
corresponding  expression  in  cylindrical  Bessel  functions,  Eq.  (2) 
is  valid  when  i  is  not  an  integer. 


2 


To  demonstrate  the  result  of  Eq.  (2),  the  method  of  Ref.  4  will 
be  followed.  If  two  functions  f(x)  and  g(x)  satisfy  the  relation 
of  Eq.  (1),  then  g(x)  is  given  in  terms  of  the  Mellin  transform  of 
f(x)  by: 


.  C+icu 

-  2nT  ic.u  x  ds/F(1's) 


(3) 


F(s)  =  ’  xS  1  f (x)dx 


(4) 


The  path  of  integration  in  Eq.  (3)  must  lie  in  the  strip  where  both 
F(s)  and  F(l-s)  are  analytic.  In  terms  of  cylindrical  Bessel  func¬ 
tions,  using  f(x)  from  Eq.  (1),  there  results: 


FU-s)  «  2  dx  J£+l/2  (x)J£+mfl/2  (x)x 
This  integral  la  a  standard  form,^  and  yields: 

. '  ( 8  +■  i);*(  i+2m+i"is) 


F(I-S>  ~sf2  — 


2&f2  s  +  m  +  s  +  l  +  m  +  s 


(5) 


(6) 


■x  m  +  1 

l  ) 


Thus,  g„  is  given  by: 
£  ,m 


g  .  _(*) 

i  ,m 


8  JL_ 
n  2ni 


ds 


2  ^  s  +  m  +  ^  s  +  £  +  "2  m  + 


i)V  - 1  ■ +  0 


x2s  :*(2S  +  i)r(x  +  \ 


m  +  I  ’  s) 


(7) 


where  s  of  Eq.  (6)  has  been  replaced  by  2s,  and  the  abscissa  of  integra¬ 
tion  lies  between  -  ^  and  0.  Regardless  of  the  parity  of  m,  the  poles 


Reference  4,  p.  214. 
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of  the  first  two  gamma  functions  in  the  numerator  are  canceled  by  the 
poles  of  f(2s  +  1)  in  the  denominator.  The  path  of  integration  may 
be  closed  by  a  large  semicircle  in  the  left  half-plane,  and  the  integral 
over  the  semicircle  tends  to  zero.  The  third  factor  in  the  numerator 
has  poles  at  s  =  -  n  -  1  +  |  m,  where  n  takes  on  all  positive  integer 
values  and  zero.  Those  poles  for  which  0  i  n  s  jm  ■  1  lie  to  the  right 
of  the  integration  path  and  do  not  contribute  to  the  Integral.  Those 
poles  for  which  j  m  s  n  i  m  •  1  have  vanishing  residues.  The  residues 
of  the  remaining  poles  may  be  evaluated,  and  the  reflection  formula 
’’(z)!'(l-z)  =  n/sin  n  z  may  be  used  to  simplify  the  expression.  Thus: 


8,  ra(x) 

l  ,m 


.  2  .  i+rof  1 

4x  (-1) 


l 


m 


IiU 

n'. 


n 


r  (2n+2-m)  xj 


2n-m 


P (r*fl-m)P(^n  +  j  -  l  -  mjP^n  +  L  +  j . 


( .s') 


0 


r(2m-n»2) 


(i«) 


2n+m 


P  (n+rof  l)P^n  +  -  ^(n  +  JL  +  m  +  j 


(9) 


Now  we  have  the  general  formula 


J  (x)J  (x) 
p,  v 


...  2rH-p,+v 

"  (-1)  P(2i>fp,+v+l)  x) 

Z,  n!  F(nHJ,+v+l)r(nHi+l)r(r>fv+l) 

0 


(  10) 


Setting  p,  =  -  A 


l  +  m  +  — ,  this  series  becomes: 


J-l-l/2(x)J£+n>fl/2('x) 


■l 


(-1)"  ['(2IH-HH-1)  (|  xj 


2n+m 


0  nl  r (nfnH-l)I'^n  +  -  fjl '^n  +■  £  +  m  +  —y 


r  no 


Reference  5,  p.  147. 
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All  the  gamma  functions  in  Eqa.  (9)  and  (11)  match  except  the  one  in 
the  numerator.  Multiply  Eq.  (11)  by  x  and  differentiate ,  and  the 
two  series  become  Identical.  Thus: 


The  identity  -jj-1/2^  an<*  return  *roTO 

cylindrical  to  spherical  Bessel  functions  now  yield  Eq.  (2).  The 
analysis  may  be  carried  through  in  the  same  manner  if  l  is  not  an 
integer.  Only  one  of  the  three  sets  of  numerator  poles  will  cancel, 
and  two  series  of  the  type  of  Eq.  (9)  result.  Both  may  be  identified 
as  Bessel  function  products,  and  combined  to  yield  the  form  of  Eq.  (2) 
in  cylindrical  Bessel  functions. 

For  large  x,  Eq.  (2)  has  the  asymptotic  form: 


**,«<*>  -  ^  -  r) 


which  includes  all  the  cases  of  Ref.  1.  For  small  x,  there  results: 


(iw-i)r  (2  +  x"*-2 

n  2- 2  r(i  +  m  +  |) 


Since  this  vanishes  for  x  ■  0,  the  calculation  of  the  potential  from 

the  phase  shift,  as  discussed  in  Ref.  1,  becomes  much  more  practical. 

The  structure  of  Eq.  (2)  makes  it  clear  that  the  inverse  function 

g.  (x)  can  always  be  written  in  terms  of  algebraic  and  trigonometric 
X,m 

functions  when  l  and  m  are  both  integers.  If  m  is  an  even  positive 

integer,  an  arbitrary  constant  may  be  added  to  g. 

»,*• 
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